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ABSTRACT 

Dark matter (DM) halos formed in CDM cosmologies seem to be characterized by a power law 
phase- space density profile. The density of the DM halos is often fitted by the NEW profile but 
a better fit is provided by the Sersic fitting formula. These relations are empirically derived from 
cosmological simulations of structure formation but have not yet been explained on a first principle 
basis. Here we solve the Jeans equation under the assumption of a spherical DM halo in dynamical 
equilibrium, that obeys a power law phase space density and either the NEW-likc or the Sersic density 
profile. We then calculate the velocity anisotropy, /3(r), analytically. Our main result is that for the 
NFW-like profile the /3 — 7 relation is not a linear one (where 7 is the logarithmic derivative of the 
density p[r]). The shape of /3(r) depends mostly on the ratio of the gravitational to kinetic energy 
within the NEW scale radius Rg- For the Sersic profile a linear /? — 7 relation is recovered, and in 
particular for the Sersic index of n w 6.0 case the linear fit of Hansen & Moore is reproduced. Our 
main result is that the phase-space density power law, the Sersic density form and the linear /? — 7 
dependence constitute a consistent set of relations which obey the spherical Jeans equation and as 
such provide the framework for the dynamical modeling of DM halos. 

Subject headings: cosmology: dark matter — galaxies: evolution — galaxies: formation — galaxies: 
halos — galaxies: kinematics and dynamics — galaxies: clusters 



1. INTRODUCTION 

In the standard cosmological model of structure forma- 
tion the luminous matter is embedded in extended dark 
matter (DM) halos. The large scale structure emerges 
out of the primordial perturbation field via gravita- 
tional instability. The model further assumes the DM 
to be made of weakly interacting particles and its dy- 
namics to be c ollisionlcss and therefore dissipationless 
(jPadmanabhani 1993 for a review). The problem of the 
dynamics of DM halos can be formulated as the classical 
gravitational TV-body problem, subject to the assump- 
tion of cosmological initial and boundary conditions. As 
such the problem can be very easily formulated, yet it 
defies any rigorous analytical treatment. The spheri- 
cal top-hat model provides the main analytical tool for 
shading light on the problem, but its scope of validity 
is rather limited (jGunn fc GottI fl972l ) . Th e model can 
be extended to accommodate shell c rossing (lGunnlll977f) 
and c osmological initial conditions ([Hoffman fc ShahamI 
Il985| ) resulting in the secondary infall model. 

In the absence of a rigorous analytical theory the study 
of the evolution DM halos relies heavily on numerical A^- 
body simulations. The advent of CPU power and im- 
proved numerical algorithm have led to a general con- 
sensus about the basic properties of DM halo, such as 
the spherically-averaged density profile, the spin and 
shape of halos. One of the pillars of the phenomenol- 
ogy of DM ha. l os is the so-called NEW density profile 
(jNavarroet al.l . ll997t ). 

It has been argued the NEW two-parameters fitting for- 
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mula provides a good approximation for DM halos found 
in a wide range of mass scales and cosmological mod- 
els, and hence can be considered universal. Subsequent 
numerical simulations have basically confirmed the func- 
tional form of the NEW profile but some controversy con- 
cerning the asymptotic slopes of the profile at small and 
large radn persists (e.g., iMoore et al.lll998l . iKlvpin et al.l 
I2OOIL iJing fc Sutdl2000[ ) . Here we refer to these general- 
izations of the NEW fitting formula as the gNEW profile. 
A second basic phenomenological finding is 
the so-called phase-sp ace density (PSD) profile 
(jTavlor fc Navarro! l200lh . These authors defined 
the PSD profile by 



Qi^) = 4: 



(2) 



where a^ is the radial velocity d ispersion. The basic find- 
ing of lTavlor fc Navarrol (|200lh is that the PSD profile 
follows a power law of the form 

Q(r)ocr-". (3) 

iTavlor fc Navarrol ()2001h found a R^ 1.875, and more re- 
cent studies find a = 1. 92 ±0.01 (iDehnen fc McLaughlinl 
|2005[) and 1.94 ± 0.01 ()Hoffman et al.ll2007[ ). Again, the 
power law behavior of Q{r) is found over a wide range of 
mass scales and different cosmological models, suggesting 
a universal property of DM halos. 

It has been recently suggested that a universal density 
slope- velocity anisotropy relat ion exists for relaxed DM 
halos ([Hansen fc Moorg l2006f ) . These author found a 
linear relation between the velocity anisotropy parameter 
f3 and the density slope 7 (< 0), 



where 



/3== -0.2(7 + 0.8), 



m = 1 - 



2cr2 



(4) 



(5) 
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and 



7 



r d p 
p d r^ 



(6) 



wfiere ag and Gcj) arc the dispersions of the two transver- 
sal velocity components. Again, the linear relation has 
been found for a variety of halos obtained in various sim- 
ulations and setups, suggesting a universal nature of the 
linear relation. 

The phenomenological gNFW density profile, the PSD 
power law and the linear /3 — 7 relation are supposed to 
provide us clues about DM halos. The universal nature of 
these relations suggests that they hold over a broad range 
of scales and models, at least within the parameter space 
of the Cold Dark Matter (CDM)-like cosmogonies. The 
equilibrium structure of collisionless self-gravitating sys- 
tems obey the Jeans equation {e.g., iBinnev fc Tremainel 
|1987(1 . which relates the velocity second moments and the 
density field. Under the assumption of a spherical sym- 
metry, the Jeans equation relates p{r), 'JT-{r) and I3{r), 
and can be rewritten so as to relate p(r), Q{r) and /3(r). 
Lacking a fundamental theory that can predict even one 
of the above relations we are motivated to study the in- 
ternal consistency of the three conditions and to find out 
whether one of these can be found to depend on the other 
two. Given the large scatter a.round the linear — 1 r ela- 
tion (|Hansen fc Moorell2006llHansen fc Stadelll2006[ ) we 
suspect that the /3 — 7 relation is the 'weakest' and less 
certain among the three relations. We therefore assume 
that the density of spherical halos follows an exact gNFW 
profile and that the PSD profile is a power law and solve 
the Jeans equation to obtain the (3 profile. This is to be 
compared with the linear (3 ~ j relation. 

It has been recently suggested that the so-called Ser- 
sic, also known as the lEinastol (Il965f). profile provides 
a better fit to the density profile (jMerritt et all [20051 : 
iPrada et al.ll2006l : iGao et all 120071 ). Given the very dif- 
ferent functional form of the Einaso profile it will be con- 
sidered here as an alternative to the NFW-like family of 
profiles and the associated /3(r) will be calculated. This 
will enable us to check the sensitivity of /3{r) to the as- 
sumed density fitting formula. 

The structure of the paper is as follows. The Jeans 
equation is solved in ij2] and the solutions of the /? profile 
are given in Sj^l and discussed in iJH 

2. THE JEANS EQUATION 

The following model of DM halo is assumed here: the 
halos are in dynamical equilibrium; it is spherically sym- 
metric; the halo density profile is given by the gNFW 
profile; and the PSD profile follows a power law. Such 
a halo should obey the Jeans equation (e.g., Binney & 
Tremaine 1987), 



d(T^.p 
dr 



MM.?,: 



-p{r) 



GM (r) 



(7) 



where M{r) is the total mass enclosed within a spherical 
shell of radius r. 

The PSD profile is expressed here in terms of the NFW 
scale radius, Rg (Eq. [12] below), the density, ps. and the 
radial velocity dispersion, a^.s-, evaluated at Rg. Namely, 



Qir) 



Ps 



R. 



(8) 



The (r/i?s)"" scaling is found to persist throughout the 
entire evolution of individu al DM halos along th e main 
branch of their merger tree ([Hoffman et al.|[2007t ). 

Rcscaling the radial coordinate to x = r/Rs the Jeans 
equation is rewritten as: 



^RsP: 



2/3 



d 
dx 



p5/3^2a/3 



_^ MMp5/3^2a/3 



GM{x) 
-P(^) ,„ ,2 -(9) 



Isolating I3{x) yields: 
/3(r) = - 



Sxp' 2a 



2/3 



GM [x) 



fT23x2"/3p2/3 R^x 



(10) 



Here the prime denotes a derivative with respect to x. 

The density is assumed to follo w a gNFW p rofile, using 
the functional form proposed bv iZhaol (|1996( ). namely 



pir) = PsPgNFW 



and 



PgNFW {x) 



2£-M 



xf^{l 



\e-A» 



(11) 



(12) 



The NFW profile corresponds to /^ = 1 and e = 3. 

Given the PSD power law and the the gNFW profile 
the structure of a DM halo is determined by three pa- 
rameters, i.e., Rs, Ps and crr,s- In particular the velocity 
anisotropy profile is given by: 

a 50^+6^) 
^^ ' 3 + 6(1 + x) ^ 

^2l+(^-M)/3^2~^/3 (1 ^ a;)2(^-A')/3 



(/i - 3) x2a/3 

2F1 (e - //, 3 - /i, 4 - /i, -x) , 



(13) 



where 2^1 (a, b, c, x) is the Gauss Hypergeometric func- 
tion and 1] = GpsR^/(J^s i^ the dimensionless con- 
stant that determines (3. Note that 77 is the square of 
the crossing-to-dynamical time ratio within Rs, hence 
it scales inversely with the (virial) ratio of kinetic-to- 
gravitational energy, T/W, within this radius. The ki- 
netic energy term, T, is calculated here taking into ac- 
count only the radial dispersion velocities, a^. 
The NFW density profile yields a /3 profile of the form 



Pnfw{x) = ~ — 



■q- 



5(l + 3x) 
6(1 + 2;) 

25/3^[ln(l + x)-y^j(l + .T) 



x4/3 



^(2q + 1)/3 



-(14) 



Hence I3{x) has a minimum around x ^ 0.1 and a max- 
imum somewhere between x ^ 1 and 10, depending on 
the value of 77. 
The Sersic density profile is given by 



where 



P{'r) = PsPSer{r/Rs) 



PSer{x) = cxp - 2n(a-^/" - 1) 



(15) 
(16) 
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Fig. 1. — The dependence of the velocity anisotropy profile /3(x = r/Rs) on the phase space density power law slope a, the inverse virial 
parameter rj = GpsR^/i^i b ^"'^'^ ^^'^ parameters of the generalized NFW density profile is presented. The left hand column displays the /3(x) 
and the right hand one the dependence of /3 on the logarithmic derivative of the density profile 7. The Hansen & Moore (2006) empirical 
fit to /3 — 7 is plotted in the right frames. A nominal model is assumed (77 = 0.16, fi = 1.0, a = 1.9 and e = 3.0) and the parameters are 
varied around that model. 



and n is the Sersic index (jPrada et al.ll2006l ). The Scrsic 
density profile yields a (3 profile of the form 



PScYix)-- 



3exp2/3 [_2ri(xi/"-l] I ^ '' 

(8"exp2/3 [-2n(a;i/" - l]:r2a/3(^^i/«)3n 
(5x^/" - a) - 2exp[2n]7ra;^/r(l + 3n) 



+6exp[2n]mTx^r]T{3n,2nx^/^))Y (17) 

where r(a;) and T{a,x) are the usual F and the incom- 
plete r functions. 

3. RESULTS 

3.1. Generalized NFW Profile 

Under the assumption of a gNFW density profile and 
a power law PSD profile the (3 profile is given by Eq. [131 
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Here we choose a base model and change the four halo 
parameters, one at the time, to check the dependence of /3 
(Fig. [T]). The base model is an N FW density profile (i.e., 
e = 3 and a /^ = 1 cusp), a = 1.9 (iDehnen &: McLaughlin! 
120051 iHoffman et al.ll200l . iRomano-Diaz et alJ (|200'7f ) 
found that T/W within i?s is roughly 2, or somewhat 
larger. Therefore, 



(18) 



-2T(< R,) ^ _3_ 

0.15. Fig. [1] shows 13 



W{< Rs) 

For the base model we choose 77 
as a function of x = r/ Rg (left column) and 7 (right 
column). The empiric al P = —0.2(7 + 0.8) relation 
([Hansen fc Moorg|2006( l is plotted for reference. 

Over the range of parameters studied here /3(.t) has 
local minimum around ~ 0.1 i?s and a local maximum 
at a fewxi?s. At the minimum the velocity dispersion 
is nearly isotropic {i.e., /? « 0. Moving outwards to 
larger radii the velocity dispersion becomes more radial, 
as P increases. However, beyond a few Rs /3 starts to 
decrease and the dispersion becomes more isotropic. It 
is interesting to note that on very small scales, below 
O.lRs, P monotonically decreases with r. Its value at 
r = depends mostly on the slope of the cusp, fi. We 
find for /i = 0.5,1.0 and 1.5 the asymptotic values of 
P for r — > are -0.21, 0.2 and 0.62, respectively. (Note 
that for fj, = 0.5 /3 reaches a minimum at r « 0.02Rs with 
P « —0.25.) The small scale behavior shows that NFW- 
like halos, obeying the phase-space density power law, 
can have steeper inner cusps by making their velocity 
dispersion more radial at the center. In other words, 
more radial orbits lead to a steeper inner density cusp. 

Eq. [131 a-iid consequently Fig. [TJ displays the depen- 
dence of /3 on ?7 — namely, an increase in 77 leads to a 
decrease in /3. The latter can be understood as following. 
Increase in rj means a decrease in T/W ratio, or equiv- 
alently an increase in (vc/cTi-)'^, where Vc is the circular 
velocity at radius r. Here two possibilities exist: (1) if 
v^ has increased at fixed a^, the system responds by in- 
creasing the tangential dispersion velocities in order to 
remain in virial equilibrium and to be supported against 
the collapse. Note that while the overall angular mo- 
mentum, J, of the system described by Eq. [7|is zero, the 
individual orbits have non-zero J and are randomly ori- 
ented which is in fact the source of the tangential velocity 
dispersion. An increase in v'^ then is translated into the 
increase in the tangential velocity dispersions in Eq. \TU\ 
and the associated decrease in p. (2) Alternatively, the 
growth in 77 can come from the decrease in a^, while w^ 
is fixed. The latter one requires that the tangential ve- 
locity dispersions stay unchanged, and, consequently, /3 
will decrease by the same token. 

3.2. The Sersic Profile 

The P{r) profile of the Sersic fit (Eq. [TT]) is presented in 
Fig. [21 The base model is taken to be ry = 0.16, a = 1.90 
and a Sersic index oi n = 6.0. The most striking feature 
of the P profile in the Sersic case is its much simpler 
functional form compared with the gNFW case. For the 
nominal a = 1.9 case the /3 — 7 relation is very close to 
linear and P grows monotonically with x over the range 
of 10~^ < X < 10. Mor eover, the Sersic index o f ti w 6 
very closely follows the [Hansen fc Moorg ([2006f ) linear 
relation. This is a surprising result. Both the gNFW and 



the Sersic profiles have been invoked as fitting parametric 
models of the density profile in the DM halos, and as 
such they do not differ substantially. Yet, their resulting 
P profiles show considerable qualitative differences. 

To gain further insight into the difference between 
the gNFW and the Sersic case we compare their cor- 
responding density (Fig. [3]) and P (Fig. \^ profiles. 
The Sersic profile has been proved to be a good fit 
for the density profile of DM halos over the range of 
10~^ < r/i?vir S, (1 ~^2)i where i?vir is the virial ra- 
dius (iPrada et al.ll2006HGao et~aI1 [20071) . The NFW fit- 
ting has been done over a comparable range. Given that 
the Sersic and the gNFW profiles are just fitting for- 
mulae and are not derived from some physical models, 
one should be cautious in using them beyond the range 
over which the fit has been performed. To prove that 
point, Figs. [31 and [H have been plotted over the range 
of 10~^ < X < 10^. The Sersic profile shows an un- 
physical drop at a; w 20 where /3 goes to minus infinity, 
namely the radial dispersion velocities vanish. This be- 
havior is shared by all the Sersic models considered here, 
of 5 ^ n ^ 8. However, DM halos do not extend that far 
and this behavior should not affect the modeling of DM 
halos. The anisotropy parameter does not converge to a 
certain limit on very small scales, x <^ 10"'^ say, and the 
monotonic dependence on the radius is not guaranteed. 

3.3. Phase Space Density of the Total Velocity 
Dispersion 

Given the radial PSD and the calculated P profile the 
PSD profile of the total velocity dispersion is easily eval- 
uated. The total PSD is given by 



Qtotir) = -3-, 



(19) 



where atot is the dispersion of the full three dimensional 
velocity field. This is related to the radial PSD Q{r) by 



Qtot(r) = Q(r)[3-2/3(r)]-3/2. 



(20) 



Assuming the nominal parameters of 77 = 0.16 and 
a = 1.90) the total PSD profile is calculated for the 
NFW and the Sersic (n = 6.0) profiles (Fig. [5]). In 
the two cases Qtot(?') is very well fitted by an r-~"*°* 
power law with atot = 1.847 and 1.837 for the Sersic 
and NFW cases, respectively. The difference in the slope 
is insignificant and we conclude that atot = 1.84. Yet, 
the fractional residual of the actual Qtot (f) from the fit- 
ted power law are as large as twenty percents for the 
NFW case but they are of the order of one percent in 
the Sersic case. The atot = 1-84 was found before by 
iDehnen fc McL aughlin' (200 51) and i s in cl ose agreement 
with the 1.82 of iFaltenbacher et all (|2007t l. 

4. DISCUSSION 

The present work has been primarily motivated by 
the qu est for understanding t he linear /9 — 7 depen- 
dence ([Hansen fc Moorg l2006( l and its relationship to 
the two so-called univers al relations, the PSD power law 
([Taylor fc Navarroll200l[ ) and the gNFW and the Sersic 
density profiles, that characterize the structure of DM 
halos. Here we have focused on whether the linear /? — 7 
relation is implied by the other two relations, and is it 
consistent with them. 
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FlG. 2. — The dependence of the /3(r) profile on a, rj and the Sersic index n is examined The structure of the figure is identical to Fig. 
[T] The base Sersic model is taken here to have 1/n = 0.16 and the lower panel shows the variation of the Sersic index. 
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Fig. 3.— The NFW (red) and Sersic (n = 6, blue) density 
profiles. The bottom panel shows the fractional difference between 
the two models. 



Fig. 4.— The NFW (red) and Sersic (n = 6, blue) (i profiles. 
In both cases r\ = 0.16 and a = 1.9 have been assumed. 
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Fig. 5. — Given the power law Q{r) profile and the calculated /3(r) the full velocity dispersion phase-space density profile, Qtot = p/u'^^^, 
is easily evaluated. Here, Qtot is evaluated for the NFW density profile, {r] = 0.16, left panel) and the Sersic model (n = 6, right panel), 
assuming a = 1.9 for both cases. The total velocity PSD closely follows a power law, Qtot oc r~"*°* , where atot = 1.84. The fractional 
residual from a power law is presented in the bottom panels. Note that in the NFW the amplitude of the residual is about 10 pcrcents and 
only about 1 percent in the Sersic case. 



We find that for no choice of parameters an NFW-hke 
density profile can yield a velocity anisotropy, /?(r), that 
is even in a rough agreement with the linear (3 — "f. This 
stands in a sharp contrast with the Sersic fitting formula 
for the DM halos density profile. For a Sersic profile 
index of n « 6.0 one recovers quite faithfully the linear 
/? — 7 relation. For that model (3{r) is a monotonically 
increasing function of the radius. For the nominal model 
we find f3{x ^ 0.1) = 0.12 and /3(x = 10) = 0.33. A 
by-product of the solution of the Jeans equation is that 
the PSD of the total velocity dispersion follows a power 
law, namely QtotC?^) oc r~"*°', with atot = 1.84, for both 
the gNFW and the Sersic density profiles. Yet, in the 
fractional deviation in NFW-like case is of the order of 
twenty pcrcents and in the Sersic case the deviation is of 
the order of one percent. 

It has been realized in recent years that the Sersic pro- 
file provides a better fit to the density profile of DM 



halos than the NFw model JMerritt et all l2005l . l2006l : 
iPrada et al.ll2006t iGao et al.lboOTt ). The present work 
substantiates and strengthens that fact and it strongly 
suggests that the the Sersic model should be used for the 
dynamical modeling of DM halos. 

The main result of the papers is that the three pillars 
of the DM halos phenomenology, namely the PSD power 
law, the Sersic density profile and the linear /3 — 7 re- 
lation, constitute a consistent set of relations that obey 
the Jeans equation. These relations provide a theoret- 
ical framework for a consistent dynamical modeling of 
DM halos. 

A very different motivation for the calculation of the 
P profile has been to provide a practical tool for mod- 
eling the mass distribution in clusters of galaxies from 
kinematic data. A powerful way of modeling the clusters 
of galaxies is based on taking moments of the velocity 
distribution of clusters' galaxies and fitting them to the 
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solutions of the Jeans equation, under the assumptions 
of the NFW density profile and /3 = const. (jLokas et alJ 
l2006HWoitak fc Lokasl[2007l ) . The presently calculated /3 
profile certainly provides a better approximation to the 
actual profile than the assumption of a constant value. 
We suggest that future analysis of the mass distribution 
of clusters will be based on the calculated profile of (3. 

We conclude the paper with a final note of caution. 
The analysis presented here applies strictly to spherical 
DM halos in virial equilibrium. Inspection of the forma- 
tion of DM halos in cosmological simulations of CDM-like 
cosmologies reveals an ongoing process where halos grow 
continuously by a slow accretion and by major mergers, 
alternating between phases of dynamical equilibrium and 
violent off-equilibrium. Therefore, one expects the DM 
halos to show some deviations from a strict virial equi- 
librium and this might limit the validity of the solutions 
of the Jeans equation. DM halos are not isolated island 
structures in an otherwise unperturbed Friedmann uni- 



verse. Halos arc experiencing an ongoing smooth, and 
occasionally no so smooth, accretion and the boundary 
between the halo and outer universe is not easily de- 
fined. The stationary Jeans equation, on the other hand 
is applicable to isolated systems. It follows that even 
in stationary, seemingly relaxed systems the Jeans equa- 
tion may not be strictly obeyed. A further complication 
arises from the deviation from sphericity of the DM halos. 
Cosmological simulations give rise to oblate and prolate 
ellipsoidal halos. This again might introduce a further 
deviation of the f3 profile of the DM halos from the so- 
lutions of the Jeans equation. These issues should be 
further investigated by means of N-body simulations. 
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